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Abstract 

Background: Malpais Spring virus (MSPV) is a mosquito-borne rhabdovirus that infects a variety of wild and feral 
ungulates in New Mexico, including horses and deer. Although, initial serologic tests and electron microscopy at 
the time of isolation nearly 25 years ago provided evidence that MSPV is a novel virus, possibly related to 
vesiculoviruses, the virus still has not been approved as a new species. 

Findings: Use of the illumina platform allowed us to obtain the complete genome of MSPV. Analysis of the 
complete 1 1019 nt genome sequence of the prototype 85-488NM strain of MSPV indicates that it encodes the five 
common rhabdovirus structural proteins (N, P, M, G and L) with alternative ORFs (> 180 nt) in the N, M and G 
genes, including a 249 nt ORF in the G gene predicted to encode a 9.26 kDa highly basic transmembrane protein. 
Although antigenically very distant, phylogenetic analysis of the L gene indicates that MSPV is most closely related 
to Jurona virus, also isolated from mosquitoes in Brazil, as well as a number of other vesiculoviruses. 

Conclusions: In sum, our analysis indicates MSPV should be classified as a member of the genus Vesiculovirus, 
family Rhabdoviridae. The complete genome sequence of MSPV will be helpful in the development of a reverse 
genetics system to study the unique aspects of this vesiculovirus in vivo and in vitro, and will assist development of 
specific diagnostic tests to study the epidemiology of MSPV infection. 
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cells, forming sizeable (3-4 mm in diameter) plaques six 
days post-inoculation [2]. Plaque reduction neutralization 
tests conducted using sera from wild mammals living in 
the area where the virus was recovered, indicated a 
relatively high prevalence (30-100%) of specific MSPV 
neutralizing antibodies among feral horses, coyotes, 
mule deer, gemsbok {Oryx gazella) and pronghorn, al- 
though the virus could not be clearly associated with 
the horse deaths [2]. 

Findings 

Methodology 
Virus source 

The prototype strain of MSPV (85-488NM) was obtained 
from the World Reference Center for Emerging Viruses 
and Arboviruses at the University of Texas Medical 
Branch. The virus stock used had been passaged four 
times by i.e. inoculation of newborn mice and three times 
in Vero cells. 
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Background 

Malpais Spring virus (MSPV) is a negative-sense RNA 
virus currently assigned to the family Rhabdoviridae in the 
order Mononegavirales [1]. MSPV was first isolated from 
pools of Aedes campestris and Psorophora signipennis 
mosquitoes collected during an investigation of unex- 
plained deaths among feral horses living near a spring in 
the Malpais Lava Flow on the White Sands Missile Range 
in New Mexico, USA [2]. Initial serologic tests and elec- 
tron microscopy provided evidence that MSPV is a novel 
rhabdovirus, possibly related to vesiculoviruses [2]. MSPV 
is lethal to newborn mice following intracerebral (i.e.) in- 
oculation and causes a cytopathic effect (CPE) in Vero 
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Next generation sequencing 

The complete genome of MSPV (85-488NM) was deter- 
mined by de novo sequencing using the Illumina platform. 

Library construction 

Viral RNA (0.05-1.7 [ig) was fragmented by incubation 
at 94° C for eight (8) minutes in 19.5 [A of fragmentation 
buffer (Illumina 15016648). First and second strand syn- 
thesis, adapter ligation and amplification of the library 
were performed using the Illumina TruSeq RNA 
Samplec Preparation kit under conditions prescribed by 
the manufacturer (Illumina). Samples were tracked using 
the "index tags" incorporated into the adapters as de- 
fined by the manufacturer. 

Sequence analysis 

Cluster formation of the library DNA templates was 
performed using the TruSeq PE Cluster Kit v3 (Illumina) 
and the Illumina cBot workstation using conditions 
recommended by the manufacturer. Paired end 50 base 
sequencing by synthesis was performed using TruSeq 
SBS kit v3 (Illumina) on an Illumina HiSeq 1000 using 
protocols defined by the manufacturer. Cluster density 
per lane was 645-980 k/mm^ and post filter reads 
ranged from 148-178 million per lane. Base call conver- 
sion to sequence reads was performed using CASAVA-1 
.8.2. Virus assembly was performed using SeqMan Laser- 
gene software (DNASTAR). In certain cases, pre-filtering 
of reads to remove host sequence enhanced the assem- 
bly process. Assembly was carried out using a fasta file 
of hamster ribosomal RNA sequences to remove host 
DNA from the assembly, thus reducing the number of 
contigs present. 

Nucleotide sequence accession number 

The GenBank accession number of MSPV prototype 
85-488NM strain is KC412247. 

Phylogenetic analyses 

The L protein sequence of MSPV was compared with 
those of 44 other rhabdoviruses downloaded from 
GenBank (members of the genera Cytorhabdovirus, 
Novirhabdovirus and Nucleorhabdovirus were excluded 
because their excessive divergence reduced phylogenetic 
resolution). All protein sequences were aligned using 
MUSCLE [3] under default settings, and ambiguously 
aligned regions were then removed using the Gblocks 
program [4]. This resulted in a final sequence alignment 
of 1136 amino acid residues. The phylogenetic relation- 
ships among these sequences were determined using the 
maximum likelihood (ML) method available in PhyML 
3.0 [5], employing the WAG+F model of amino acid 
substitution and subtree pruning and regrafting (SPR) 
branch-swapping. The robustness of each node was 



determined using 1,000 bootstrap replicates utilizing 
nearest neighbor interchange (NNI) branch-swapping. 

Antigens and immune reagents 

Antigens used in complement fixation (CF) tests and for 
immunizing animals were prepared from infected new- 
born mouse brains by the sucrose/acetone extraction 
method [6]. Hyperimmune mouse ascitic fluids were 
prepared against each of the 21 rhabdoviruses listed in 
Table 1. The immunization schedule consisted of four 
intraperitoneal injections given at weekly intervals. Im- 
munogens consisted of 10% crude suspensions of ho- 
mogenized infected mouse brain mixed with equal 
volumes of complete Freund's adjuvant just prior to in- 
oculation. Sarcoma 180 cells were also given intraperito- 
neally after the final immunization in order to induce 
ascites formation. All animal experiments were carried 
out under an animal protocol approved by the University 
of Texas Medical Branch lACUC committee. 



Table 1 Antigenic relationships of MSPV with confirmed 
and tentative members of the genus Vesiculovirus 



Genus 


Viral antigen 


MSPV 
antibody* 


Homologous 
antibody 




Carajas (CSV) 


0 


128 




Chandipura (CHPV) 


0 


256 




Cocal (COCV) 


0 


256 




Isfahan (ISFV) 


0 


64 


Vesiculovirus 


Maraba (MARV) 


0 


512 




Piry (PIRV) 


0 


64 




Moreton (MORV) 


0 


256 




VSV Indiana (VSIV) 


0 


256 




VSV New Jersey 
(VSNJV) 


0 


32 




Boteke (BTKV) 


0 


256 




Farmington (FARV) 


0 


256 




Gray Lodge (GLOV) 


0 


8 




Jurona (JURV) 


0 


1024 


Tentative 


KlarTiath (KLAV) 


0 


256 


Vesiculovirus 


Kwatta (KWAV) 


0 


64 




La Joya (LJV) 


0 


512 




Perinet (PERV) 


0 


512 




Porton (PORV) 


0 


256 




Radi (RADV) 


0 


512 




Jug Bogdanovac 
(YBV) 


0 


256 




Malpais Spring 
(MSPV) 


512 





* Antibody titers expressed as reciprocal of highest positive antibody dilution. 
0 = <8. 
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Serologic testing 

Complement fixation tests were conducted according to 
a microtechnique described previously [7], using 2 full 
units of guinea-pig complement. Titers were recorded as 
the highest dilutions giving 3+ or 4+ fixation of comple- 
ment on a scale of 0-4+. 

Results 

The complete genome of MSPV (85-488NM) showed that 
the 11,019 nt MSPV genome exhibits typical rhabdovirus 
organization and consists of 3'- and 5'- terminal non- 
coding regions (approximately 63 nt and 97 nt, respec- 
tively) and 5 coding regions (genes), each bounded by pu- 
tative transcriptional regulatory sequences (Figure lA). 
Similar to those of most other animal rhabdoviruses, the 
transcription initiation (AU[U/C]GUC) signal and tran- 
scription termination/polyadenylation ( [G /U] UAC [U] 7) 
signal are relatively conserved (Figure IB). Gene 1 (N), en- 
codes a 426-aa protein (47.4 kDa); gene 2 (P), encodes a 
274-aa acidic protein (29.8 kDa; pi 4.75); gene 3 (M), en- 
codes a 226-aa protein (25.6 kDa); gene 4 (G), encodes a 
522-aa protein, which is a class 1 transmembrane glyco- 
protein; and gene 5 (L), encodes a 2096-aa protein (238.7 
kDa), which the RNA-dependent RNA polymerase (RdRp). 
Alternative ORFs (< 180 nt) also occur in the N gene (198 
nt), M gene (189 nt) and G gene (249 nt) encoding po- 
tential polypeptides for which BlastP searches revealed 
no significant sequence relationships to known proteins 



(Figure lA). The putative proteins encoded in the N and 
M genes are unremarkable in structure. The putative 
83-aa product encoded in the G gene is predicted (Phobius 
and HMMTOP) to be a highly basic transmembrane pro- 
tein (9.26 kDa; pi 11.36) with a 53-aa N-terminal cytoplas- 
mic domain, a 24-aa helical transmembrane domain and a 
very short (6-aa) C-terminal ectodomain (Figure IC). 
ORFs encoding similar small hydrophobic proteins of un- 
known function have been identified in the genomes of 
other rhabdoviruses [8]. The MSPV G protein contains 12 
cysteine residues in the ectodomain that form six disul- 
phide bridges that are common to all vesiculovirus G pro- 
teins and three predicted N-glycosylation sites, two of 
which align with those of the vesicular stomatitis Indiana 
virus (VSIV) G protein [9]. 

Serologic studies conducted at the time of MSPV iso- 
lation indicated a weak antigenic relationship with an- 
other probable vesiculovirus Jurona virus (JURV) by 
indirect immunofluorescence, but none with VSIV, ves- 
icular stomatitis New Jersey, Cocal, Piry or Chandipura 
viruses (CHPV) [2]. We could not demonstrate an anti- 
genic relationship in CF tests between MSPV and an ex- 
panded set of confirmed and tentative vesiculovirus 
species, including JURV (Table 1). 

Amino acid sequence alignments (MUSCLE) using 
available rhabdovirus N, G and L proteins indicated that 
MSPV shares highest sequence identity to vesiculoviruses. 
Phylogenetic analysis based on an alignment of 1136 
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Figure 1 Genome organization of MSPV. A. Schematic representation of tine genome organization of MSPV. Blocl< arrows indicate tlie location 
of long ORFs, including overlapping ORFs in the N, M and G genes (light shading). B. Sequences of MSPV gene junctions illustrating potential 
transcription initiation sequences (shaded), conserved transcription termination/polyadenylation sequences (solid underlined) and initiation 
codons (dotted underline). C. Schematic representation of the predicted membrane topology of a small, basic transmembrane protein encoded 
in the alternative ORF in the G gene. 
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Figure 2 Phylogenetic tree of L protein sequences from 45 rhiabodoviruses. The position of MSPV is sliaded and bootstrap support values 
(>70%) are sliown for l<ey nodes. All horizonal branch lengths are drawn to a scale of amino acid substitutions per site, and the tree is rooted in 
the position observed in a broader analysis of the Rhobdoviridoe (i.e. including members of the Cytorhobdovirus, Novirhabdovirus and 
Nucleorhabdovirus; not shown). (*) Genera not formally approved by the International Committee on Taxonomy of Viruses i\CT\/). 



amino acid residues of the L protein sequence of 45 rhab- 
doviruses indicated that MSPV is most closely related 
(although with weak bootstrap support) to JURV, which 
has been isolated from Haemagogus spp. mosquitoes in 
the state of Para, Brazil in 1962 [10,11], as well as to 
Isfahan, CHPV and Perinet viruses (Figure 2). These five 
viruses form a well-supported (100% bootstrap support) 
phylogenetic cluster within the genus Vesiculovirus, 

Discussion 

The genome organization and phylogenetic relationships 
indicate that Malpais Spring virus should be classified as 
a new species in the genus Vesiculovirus, This is also 



consistent with the previous report of a distant antigenic 
relationship with JURV, which also falls phylogenetically 
within the vesiculovirus cluster. Evidence of alternative 
ORFs in the MSPV N, M and G genes is unique 
amongst vesiculovirus genomes sequenced to date but 
overlapping ORFs appear to occur commonly in rhabdo- 
viruses [8], including VSIV in which an alternative ORF 
in the P gene encodes the 55-aa and 65-aa C and C pro- 
teins which are known to be expressed in infected cells 
from alternative initiation codons [12,13]. The complete 
genome sequence of MSPV will be helpful in the deve- 
lopment of a reverse genetics system to study the unique 
aspects of this vesiculovirus in vivo and in vitro, and will 
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assist development of specific diagnostic tests to study 
the epidemiology of MSPV infection. 
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